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The stochastic Eden model of charged particles aggregation in two-dimensional systems is pre- 
sented. This model is governed by two parameters: screening length of electrostatic interaction, 
A, and short range attraction energy, E. Different patterns of finite and infinite aggregates are 
observed. They are of following types of morphologies: linear or linear with bending, warm-like, 
DBM (dense- branching morphology), DBM with nucleus, and compact Eden- like. The transition 
between the different modes of growth is studied and phase diagram of the growth structures is 
obtained in co-ordinates. The detailed aggregate structure analysis, including analysis of their 
fractal properties, is presented. The scheme of the internal inhomogeneous structure of aggregates 
is proposed. 

PACS: 02.70.Lq, 68.70.+W, 75.40M 



I. INTRODUCTION 



There are a lot of computer models aimed at simulation of the natural patterns occurring in colloidal aggregation, 
electrochemical deposition, dielectric breakdown, bacterial colony growth, viscous fingering, spreading, etc. Such 
models allow to explore the process of pattern formation in real physical systems and they are based mostly on the 
Eden model |^ and the model of diffusion limited aggregation (DLA) with due account of particles drift, convection 
flow, diffusion, external field influence and specific interactions between the particles. The random successive growth 
model two-scale drift-diffusion model diffusion- migration and convection- migration models are good 
examples of such models. These models describe the most important pattern morphologies observed in various non- 
equilibrium systems, such as DLA-like, dendritic, needle, treelike, dense-branching, compact, stingy, spiral and chiral 
structures 

The relationship between the structure and the growing pattern morphologies, growth mechanisms and peculiarities 
of interaction between the elements of growing structure seems to be the main problem in this field. The local 
interaction is introduced into the screened growth model ||], in which the growth probability is determined by 
multiplicative screening contributions of the occupied sites nearest to the site of growth. In this model, the fractal 
dimensionality of the cluster Df is determined by the form of the screening function and may be varied in the wide 
range. Ausloos et al. have studied the magnetic Eden model Q and magnetic diffusion-limited aggregation model 
[ |TT[ for the Izing-like nearest-neighbour interaction between the spins. Each of these models leads to the formation of 
different types of aggregates for different values of parameters controlling the growth processes. However, inclusion 
of the long-range particle interactions, for example, of dipolar type, requires considerable CPU-resources. At present, 
such simulations are done only for the systems of moderate size However, this question is very important for the 
interpretation of the experimental data on magnetic particles aggregation [|l^ . 

The study of the models of charged particles aggregation is also of a great interest. It is especially true in relation 
with the study of particle surface charge effects on the mechanism of colloidal aggregation | 

in electrochemical deposition experiments [|l^, pattern formation mechanism in discharge systems [|16[ and possible 
fractal structure of ball lightning [|l^ . The presence of the long-range interactions can produce the qualitatively new 
effects, not usual for the systems with the short-range interactions. Thus, systems with the spatial charge distribution 
are characterised by the limited stability. It is well known, that the electrostatic interactions cause the loss of stability 
of charged droplet [l8|, and may cause instabilities of charged polymer systems aerosols , etc. For the non- 
equilibrium growth condition, one may expect the occurrence of such phenomena as large-scale fluctuations, lacunarity 
or growth of finite size aggregates. 
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In this work, we propose the stochastic Eden-hke lattice model of charged particles aggregation. Particles overcome 
the electrostatic repulsive barrier created by the aggregate and stick to it due to existence of a short-range attraction 
force. We observe the formation of aggregates with different morphologies, such as: linear or near-linear, linear with 
bending, warm-like, DBM or DBM with nucleus, and compact Eden-like. The transition between the different modes 
of growth is governed by the two model parameters screening length. A, and attraction energy, E. 

The structure of the paper is as follows. In Section || we describe the model. Section |l| reports the results 
of simulation and their discussion, including examples of the clusters, phase diagram, analysis of the radial cluster 
structures and their fractal properties. Section |v| contains our conclusions. Appendix ^ presents some simple 
electrostatic estimations for the model. 



II. MODEL 

The model is based on the standard two-dimensional stochastic Eden model We assume that all the particles 
bear the same electric charge q. A particle located at point j in the proximity of the aggregate interacts with all the 
particles of aggregate through the long-range Coulomb repulsion, > 0, and short-range 'glue' attraction, E. We 
also take into account the possibility of screening of electrostatic repulsion, and introduce the screening length A, in 
order to switch off the influence of the aggregate particles located further then at the distance A from point j. The 
energy of electrostatic repulsion may be defined as: 

where r^j is a distance between the j-point and z-particle of the aggregate, aij = 1 for r^j < A and a^j = in any 
other cases. The cluster growth takes place on the simple square lattice. The lattice step and the particle charge are 
assumed to be equal to 1. So, in this model, the cluster growth is governed by two parameters: the screening length, 
A, and the attraction energy, E and we may consider the values of A, U and E as dimensionless. 

A cluster grows according to the following scheme: 

a) the seed particle of the cluster is placed at the lattice point x = y = 0; 

b) the cluster perimeter site, j, is chosen randomly (the perimeter of the cluster is a set of unoccupied lattice sites 
nearest to the cluster); the value of U for the j-site is calculated and compared with the value; 

c) if U < E^ then the j-site becomes the site of a cluster, and determination of a new cluster perimeter site takes 
place; otherwise, if U > E^ the j-site is declared as unoccupied forever; 

d) then the step b) is reiterated. 

These steps run until the cluster reaches either the lattice size, or the desirable size, or until the cluster growth 
termination. The last case may be observed when the following condition becomes true for all the sites of cluster 
perimeter: 

U> E. (2) 

The lattice size varied from 400x400 to 5000x5000, the maximal number of particles in the cluster did not exceed 
310^ The results were averaged over the 10 ^ 30 different initialisations. 



III. RESULTS AND DISCUSSIONS 



A. Phase diagram 

Fig. |l| presents examples of clusters that appear at different values of the screening length. A, and attraction energy, 
E. All patterns were obtained on the basis of cluster images on 400x400 lattice. Small circle in the centre of each 
pattern denotes the seed particle. The plane of A and E values is approximately divided into several regions that 
correspond to the observed cluster types. These regions are schematically presented on the phase diagram in Fig. ||. 
The numbers of patterns of the Fig. ^ correspond to the numbers in black squares put on the phase diagram. We 
discern the following regions of the phase diagram: 

Lf is the region of linear or near-linear finite cluster growth; 

Lboo is the region the infinite linear clusters or linear clusters with bending; 

Wf is the region of finite warm-like clusters; this type of structure is represented by the patterns 1, 5, 9, 13, 17 of 
Fig. ^ and similar structures are reported to occur for the random successive growth model B; 
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Woo is the region of infinite warm-like clusters; this structure type examples are presented by the patterns 2, 6, 10, 
14 of Fig. I; 

DB oo is the region of lacunary clusters, clusters with developed branch structure and clusters of DBM-type H] 
structures growth; in this region, clusters may possess a central dense part or a nucleus (see, patterns 3, 4, 7, 8, 11, 
12, 15, 16, 18, 19, 20 of Fig. 0); 

Coo is the region of compact clusters or compact clusters with cavities growth; in the limit E ^ oo these clusters 
are the clusters of standard Eden model. 

The case of A = oo (see 21-24, Fig. 0) is shown above the main diagram (Fig. ^). There are no regions of infinite 
cluster growth and no compact cluster growth for A = oo (except for the limiting case of £^ = oo, when the growth of 
Eden-like clusters takes place). 

Thus, on increasing of E at fixed A value, we observe the following sequences of regions 

Lf^ Lboo^ Wf ^ Woo^ DBoo^ Coo (3) 

when A is finite, and 

Lf ^ Lbf ^ Wf ^ DBf (4) 

when A=oo. 

The main distinctive feature of the phase diagram in Fig. || is the existence of regions with finite and infinite 
cluster growth. It appears merely due to the long-range character of Coulomb interactions. When a cluster grows, 
its total charge increases and it results in increase of the electrostatic repulsion energy at the points of new particle 
attachment, which are located on the cluster perimeter. In the certain cases, when the cluster sufficiently large grows, 
the condition of cluster growth termination (||) realises for all the points of cluster perimeter. We characterise the 
maximal size of a cluster by its maximal radius i^max, which is simply the distance from the central seed particle to 
the most remote particle of the cluster. The full manifestation of the long-range character of Coulomb interactions 
occurs only in the cases when the value of A exceeds the spatial dimension of the cluster, i.e. when A > i^max- This 
condition is always true when A = oo, which is the cause of only finite cluster growth. 

How it is possible to discern the different regions of the phase diagram presented in Fig. ||? Some conclusions about 
localisation of different regions can be made on the base of the visual analysis of patterns similar to that presented on 
Fig. [l[ For fixed A, the increase of E leads to the growth of different type clusters. Thus, in Lf region, the length of 
cluster increases with E and in Lb region the number of cluster bendings increases and the length of the linear cluster 
part decreases with the increase of E. However, the visual analysis does not allow us to make an exact estimation of 
region localisations. This method is effective only for determination of the sharp boundary between the region DBoo 
and region Coo- Note that this boundary can also be localised by means of simple electrostatic estimations from 
which we obtained for the transition energy Ec between region DBoo and region Coo the following relation: Ec = ttA 



(See Appendix, Eq. (|A5|)). 

The boundaries of Lf , Lboo, Wf and Woo regions are not defined clearly and the transition from one phase to 
another is not sharp and appears to be rather smeared. Due to the finite-size effects near the boundaries of regions, 
we observe strong ffuctuations of cluster sizes and morphologies. The width of such transition zones between the 
different regions may vary and is difficult to determine, however, in our simulations it does not exceed AE ^ 0.2 for 
fixed A. On moving from region Lf to region Loo at fixed A value, the closer we are to region Loo, the higher is the 
probability of an infinite linear cluster growth. Correspondingly, the probability of finite linear cluster growth is lower. 
For estimation purposes, the boundary line may be defined as a line of approximately equal probabilities of infinite 
and finite clusters formation. This is true also for the cases of Lboo— Wf and Wf — Woo boundaries. Unfortunately, 
such approach is not acceptable for boundary localisation between Woo and DBoo regions. 

To estimate the region boundaries more precisely, we have carried out the configuration analysis of cluster structures 
and studied the behaviour of the maximal cluster radius R ^ . 



B. Configuration analysis and maximal cluster radius 

Transitions between the different regions of phase diagram may be also analysed with the help of configuration 
analysis. The idea is to find some correlation between the pattern type and the probability of certain nearest neighbour 
configuration appearance in that pattern. Any cluster on a simple square lattice consists of the particles that have 1, 
2, 3 and 4 nearest neighbours. For the purposes of configuration analysis, we have discerned four types of configuration 
presented in Fig. ||. 

It is easy to see that the configurations of type (1) correspond to the points of cluster growth termination, the 
linear clusters consist mostly of type (2) linear configurations, the clusters with bending and warm-like clusters are 
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the mixtures of type (2) and type (3) configurations of particles and the particles of compact and branched clusters 
form mostly the configurations of type (4). 

The fractions of particles in (i) type configurations, P^, were calculated as the ratio of the number of particles 
having one neighbour in configuration (1), two neighbours in configuration (2), two neighbours in configuration (3) 
and three or four neighbours in configuration (4) respectively, to the total number of cluster particles. 

To localise more precisely the boundary between the above mentioned finite and infinite cluster regions, we have 
also studied the dependency of maximal cluster radius i?^^^ versus attraction energy E for fixed values of screening 
length A. The results were averaged over 10 ^ 30 different initialisations for each chosen pair of values A, E and 
the maximal lattice size used this work was 5000x5000. A cluster reaching the lattice boundaries was considered as 
infinite. Such clusters were observed only for A < oo, while, as it is mentioned above, for A = oo all clusters are only 
finite, they stop to grow far beyond the limits of our system. 

Figures |^|| show the fractions of particles with relevant nearest neighbour configurations. Pi (a) and the maximal 
cluster radius, Rmax versus E (b) for A = 30, oo, respectively. These figures correspond to the continuus move along 
the E axis across the phase diagram presented at Fig. || at fixed values of A. 

Note that at finite A the transition between the finite and infinite growth regions in the sequence (H) is obvious at 
^max versus E curve (Fig. ^). This is the way to fix the transition boundaries between the Lf — Lboo, Lb^— Wf , or 
Wf — Wcx) regions quite reliably. 

In the region of linear or near-linear finite cluster growth (Lf ) the considerable increase of P2, is observed with E 
increase with relatively weak Pi, P3 and P4, vs. E dependences. Such P2 increase is accompanied by the increase 
of R^^^. The Lj— Lboo boundaries (for A = 30, Fig. ||) or Lf— Lbf boundaries (for A = 00, Fig. ||) at £^ = E}j{X) 
approximately correspond to the region of maximal P2 values. 

We can estimate the R^^^{E) dependence in the Lf region, as well as the E = £^5 (A) boundary between the Lf and 
Lboo regions (at finite A) proceeding from the following electrostatic considerations. The mainly linear clusters grow 
at small values of E in the Lf region. In this region the particles tend to join the chain cluster at the tips where the 
repulsive energy U is the lowest. The particles forming a linear chain of length L , repulse with the energy U the 
newcomer at the edge of the chain: 

U (X ln(A) (5) 

when A(< L) is finite, and 

U (X ln(L) (6) 

when A = 00. 

The growth of the linear chain terminates, when the cluster growth termination condition (H) is realised, so that 
we have for the maximal radius at £^ < Eb{X): 

i^max cx exp(£;) (7) 

where, the (A) dependency may be approximated for finite A(A < 200) by 

^5(A)^1.21n(A). (8) 

which corresponds to the boundary line between the Lj and Lboo regions at the diagram of Fig. H and Eb{X = 
00) = 8.0 is the boundary line between the Lf and Lbf regions. 

Note, that at any A value the Eq.(|^) approximates quite well the behaviour of R^^^{E) in the interval of £^ < £^5 (A) 
see Figs. |b-|b). 

The transition to the region W/, where the growth of finite warm-like clusters takes place (both for A = 30 and 
for A = 00), is accompanied by the sharp decrease of R^^^ with the simultaneous increase of the fraction of particles 
in type (3) configurations, P3. We approximately localise the Woo— DBoo (for A = 30) or Wf — DBf (for A = 00) 
boundaries for the transitions between the warm-like and DBM structures at those values of E^ which correspond to 
approximately the same fractions of particles in configurations of (2), (3) and (4) types (P2 ^ Ps ^ P4 ^ 0.2 ^ 0.3). 
Note, that in all the cases the maximum point of Pi is localised near these values of E. When E is sufficiently large, 
the type (4) configurations prevail in the transition zone between the Woo and DBoo regions. 

C. Radial structure of cluster 

The mechanisms of the internal cluster structure formation, as well as the matters dealing with its morphology may 
be clarified by the analysis of certain cluster properties changes of in the radial direction. For this purpose, we have 
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analysed the local density profiles p{r) and the E/U{r) ratio that characterise both the cluster growth stability and 
the possibility of cluster growth termination if the condition (H) is realised. 
We define the local density profile as 

p(r) = N{r)/N,{r), (9) 

where N{r) is the number of cluster particles inside a ring of radius r and one lattice unit width, centered at a seed 
particle; Nc{r) is the similar number of densely packed cluster particles. Such definition implies that the denser is the 
cluster, the closer is its density to 1. 

Consequently, according to such definition, the value of p corresponds to the relative density profile reduced to 
a density profile of densely packed cluster. This definition allowed us to suppress the long-range effects of magic 
numbers, that are clearly present in the density profiles of densely packed clusters on regular lattices 

Figures ^js] present the profiles of local density, p (a) and E/U (b) for clusters obtained at A = 30, 100 and oo, 
respectively. These dependencies were obtained at E = 10, 20, 100 for each A and averaged over 10 clusters grown 
on 200x200 lattice (for A = 30, 100) and 150x150 lattice (for A = oo). Thus, we have analysed the structure of clusters 
in the following regions: Woo,DB^ and Coo (for A = 30), Woo and DBoo (for A = 100), as well as the regions Wf 
and DBf (for A = oo ) (see the phase diagram in Fig. ||). 

In all the cases we have observed a rather sharp and self-correlated initial decrease of p and E/U values with r 
increase. For A = 30 and A = 100 these dependencies have minimum at r ^ A/2 and maximum at r ~ A. Analysis of 
the p{r) and E/U(r) profiles allows us to make some general conclusions dealing with the character of internal cluster 
structure inhomogeneity. In the most general case, four different zones of spatial cluster structure may be discerned 
along the radial direction of cluster structure. (Fig.^: 

zone A or internal compact nucleus; this nucleus zone has the density p ~ 1.0 and its radius may be estimated as 



Rn ^ E/4 (See Appendix, Eq.(Q); 

zone B or near- nucleus boundary shell; in this shell, the density drops approximately as p oc 1/r or all the cases 
observed, and this shell is located in the region of Rn < r < A/2; the existence of this shell is connected with the 
growth of near-linear branches from the dense nucleus; 

zone C, or screening- induced compression shell; this shell is located in the region of A/2 < r < 1.2A and may be 
present only for finite A's; the density maximum here is observed at r ~ A; this shell exists due to the cluster structure 
compression effects, resulting from certain screening limitation of electrostatic repulsion of particles from the most 
dense part of cluster (its nucleus) at distances within the above-mentioned range; correspondingly, the maximum of 
the E/U{r) curve is observed in this zone; 

zone D, or external zone; this zone is characterised by the gradual density decrease and, in some cases, the fractal 
behavio ur of density is observed: p oc l/r^~^^, where Df is the fractal dimensionality of cluster in this zone (See 
Section pIDl ). 



For the case of A = 30 (Fig. |6|), the energy, E^ increase from 10 to 100 corresponds to the gradual transition 
Woo DBoo Coo on the phase diagram. At £^ = 10 (the lowest curve), we are at the border of Woo DBoo 
(for example, see Fig. I, cluster 11). The zone A is absent in this case, since the nucleus radius is very small 
(Rn ^ E/4 = 2.5). At 20 < £^ < 80, we are in DBoo region (for example, see Fig. cluster 12); here, each energy, 
E, increase by 10 results in a nucleus size increase, approximately, by 2.5. In these cases, the profiles of p{r) and 
E/U{r) reveal somewhat nontrivial behaviour, and we can find the zones A ^ D in the cluster structure. At £^ > 80, 
we are almost near the border of Coo region. We observe that p{r) and E/U{r) profiles elevate with E increase; these 
profiles decrease with r without passing through maxima and minima. The dashed lines on Fig. ^ correspond to the 
case of £^ = 130, i.e., when we are far in the Coo region, where the clusters are densely packed or compact. In this 
last case, p{r) ^ 1.0 is observed and the different zones of cluster structure may be discerned only by the character of 
E/U{r) decrease. 

The similar behaviour of p{r) and E/U{r) profiles is observed also in the case of A = 100 (Fig. I). However, 
some aspects of the cluster structure are more explicit in this case. Here, we have the more prominent maxima and 
minima of p(r) and E/U(r) profiles, and we observe no general elevation of E/U{r) curves with E increase. It may 
be explained by the fact that such behaviour is pronounced only at sufficiently large values of E in the proximity 
of compact growth region Coo- We had not reached Coo region in our simulations: for A = 100, it starts only at 
E > 300. 

At A = oo, the profiles of p{r) and E/U{r) (Fig. |) do not pass through maxima and minima, but simply drop 
almost like 1/r outside the nucleus zone, A. Figures l|, 21-24 contain good examples of such clusters. The dashed 
lines (Fig. H) correspond to the cases when E = 200(1) and 300(2). In these cases, the clusters have the large compact 
nuclei with the radii of Rn ^ ^0 — 75. 

We have also analysed the behaviour of the E/U{r) profiles for sufficiently small values of E in such regions as Lf, 
Lboo and Wf for A = 30 (Fig. |lo| a) and in such regions as Lf , Lbf and Wf for A = ooO (Fig. p!Q|b). Each of these 
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E/U{r) profiles was obtained through data averaging over the 10 different clusters. In the regions of finite cluster 
growth, the E/U{r) decreases quickly with r upraise. At certain critical value of r = R^^^^ we observe E/U{r) = 1 
at all the points of cluster perimeter and just at this stage the cluster stops to grow. In transition zone between the 
finite and infinite growth regions, E/U{r) takes a certain stationary value greater than 1. 

Remarkable is the case of A = 30 and E = 4.5 (Fig. p!o|a), when we get into the narrow Lboo region on the phase 
diagram of Fig. |2[ We see that after the fast initial decrease of E/U{r)^ its abrupt uprise is obrerved at the certain 
r ~ A, and then the E/U(r) value of does not change with r. This condition corresponds to formation of the infinite 
linear clusters or linear clusters with bending. 



D. Fractal properties of clusters 



As we have shown in the previous Section [III C| , in the general case, the clusters obtained within our model do 
not reveal any homogeneous properties in their radial direction. Therefore, the fractal approach is not applicable 
for analysis of the structure of whole cluster. Instead, we have analysed the fractal properties of cluster branches in 
the zone D, or its external zone, where we observed the gradual local density decrease, which may be approximately 
described by the relation p oc l/r^~^^, where Df is the fractal dimensionality of cluster in this local zone. 

In order to estimate the Df value in this zone, we apply the sandbox method with shells of increasing outer radius, 
excluding the central cluster part of the radius R ^ 1.2 A. We analyse the dependency of the number of particles in 
the shell, versus its area, using the following relation: 

N (X S^f^^, (10) 

The results are averaged over 30 clusters for each value of {E^ A) pair. 

Figure |ll| presents the results of Df calculations for A = 10,30,50, 100 versus E/X. On E increase we observed 
the transition from linear {Df = 1) to compact {Df = 2) structures. The detailed analysis shows that this transition 
occurs when the system is in the region where the growth of finite (W/ ) or of infinite (Woo) warm-like clusters takes 
place. Although the fractal dimensionality Df can be slightly lower than 2 in the DB^o region, however, it approaches 
to this value as the system approaches the Coo region, where the clusters are dense and their fractal dimensionality 
Df is strictly equal to 2. For the case when A = oo the branches of DBf clusters do not display uniform fractal 
properties, but their fractal dimensionality tends to 1 when we move in radial direction from the centre. So, in this 
case we have Df = 1 and this behaviour is easy to understand in view of the tendency of electrostatic repulsion to 
linearise the branches of the cluster. 

We think that the smooth character of Df transition on passing through the regions W/ or Woo results from the 
stochastic nature of the model we used in our simulations. Quite different behaviour of fractal dimensionality can be 
observed for the similar, yet deterministic in nature model of aggregation [^. In this case we observed the sharp 
transition between Df = 1 for finite clusters and Df = 2 for infinite clusters. 



IV. CONCLUSIONS 



We have introduced the two-parametric electrostatic Eden model that simulates the aggregation of charged particles. 
Cluster growth is controlled by the screening length of electrostatic interactions A, and the energy of particle adhesion 
to a cluster E. The phase diagram in A versus E co-ordinates reveals the regions of finite (f) and of infinite (oo) 
cluster growth. The model gives a variety of patterns with different morphology, including fractals. We do discern 
the following morphologic types: linear, linear with bending, warm-like, dense-branching and compact Eden-like. 
The structure compactness increases with the attraction energy E increase, while the strengthening of electrostatic 
repulsion on A increase results in linear morphologies predominance. The cluster structure is not homogeneous in 
the radial direction and in the general case four zones with different morphological properties may be differentiated: 
internal compact nucleus zone, near-nucleus boundary shell, screening induced compression shell, and external shell. 
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APPENDIX A: SIMPLE ELECTROSTATIC ESTIMATIONS FOR THE MODEL 



We approximate a compact two-dimensional Eden cluster by uniformly charged disk of radius R. We calculate the 
energy of electrostatic interaction between the charge and compact disk-like cluster of charged particles U by using 
numerical integration instead of summation in Eq. ( |l|) and taking into account the homogeneous charge density 
distribution in the disk-like cluster with p{r) = 1. 

The solid line in Fig. |l^ presents the dependency of the reduced electrostatic energy, [/* = U/R, against the 
reduced distance between the charge and the cluster centre, r* = r/R, for A = oo. It is obvious that for the particle 
located at the cluster centre, when r = 0, the energy is equal to Uq = 27tR. At large distances from the cluster centre, 
when r* ^ oo, we may neglect the spatial charge distribution inside the cluster and treat the whole charged cluster 
as a charge Q = ttR'^ concentrated at the centre of such cluster. In this case: 

[/*(r* -^oo) = U^= 7r/r*. (Al) 

The dependency of on r* is presented on Fig.O by the dashed line. We see, that, at least, for r* > 1.5, the 



Eq.(Al) approximates the /7*(r*) dependency rather precisely. 

If the particle is located at the circle-like cluster perimeter, where r = (r* = 1), the interaction energy expressions 
are easy to obtain in the analytical form for the arbitrary A value: 



[/; = 4(1 + A* arccos(A*) - Vl - A*^) for A* < 1 (A2a) 
U; = 4 for A* > 1 (A2b) 

Here, A* = \/{2R). 

The condition A* > 1 from Eq. ( |A2b|) is equivalent to the condition A > 2R. In this case, the Up value is defined by 
the interactions with all the cluster charges, without any screening limitations. Thus, for all the R < A/2, we obtain 



from Eq. (|A2b| ): 

Up = AR. (A3) 

The stable growth of compact aggregates terminates when E < Up] therefore, in the specific case of A = oo we may 
estimate the radius of the compact nucleus as 

Rn ^ E/A, for A = oo (A4) 

At R > Rn, the loss of compact cluster growth stability takes place and cluster ramification begins. 



It is possible to estimate the Rn values for finite screening lengths, A's, by substituting Up = E into Eq.( |A2a| ). 



However, in this case the situation of infinite compact cluster growth with Rn = oo is possible for sufficiently large 



values of E > Ec. For the purposes of Ec estimation, we consider Eq.(A2a) when A* ^ (or R ^ oo when A is a 



finite) and obtain U^ = 27rA* or Up = ttA. So, the relation for Ec{X) may be written as follows: 

Ec{X) =7rA. (A5) 



The Eq.(A5) allows us to estimate the position of the boundary between the Coo and DBoo regions on the phase 
diagram presented in Fig.^. 
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FIG. 1. Examples of cluster patterns at different values of the attraction energy E and screening length A. Numbers of 
patterns correspond to the numbered black squares on the phase diagram on Fig.g. All patterns are displayed on the same 
scale 500x500. 

FIG. 2. Phase diagram for the cluster growth in co-ordinates of screening length A versus attraction energy E. See the text 
for details. 

FIG. 3. The four types of the nearest neighbour configurations of cluster growing on a square lattice discerned for the 
purposes of configuration analysis. 

FIG. 4. The fractions of different types of particles configurations, Pi, i = 1 ^ 4 (a), and maximal cluster radius, i?max,(b) 
versus attraction energy, ^, at A = 30. The value of Eb determines the boundary between the Lf and Lboo regions and 
may be approximated by Eq.(^. The value of Ec determines the boundary between the DBoo and Coo regions and may be 
approximated by Eq.( |A5| ). The dashed line corresponds to the Rma^{E) curve in the interval of E < Eb obtained through rms 
approximation using Eq.(^). 

FIG. 5. The fractions of different types of particles configurations, P^, i = 1 ^ 4 (a), and maximal cluster radius, i?max ,(b) 
versus attraction energy, ^, at A = oo. The value of Eb{^ 8) determines the boundary between the Lf and Lbf regions. The 
dashed line corresponds to the Rma^{E) curve in the interval oi E < Eb obtained through rms approximation using Eq.(0). 

FIG. 6. The radial profiles of local cluster density, p(r), and the energy ratios, E/U{r). The results correspond to the 
different values of attraction energy E = 10,20,30,40,50,60,70,80,90,100 (solid lines), E = 130 (dashed line), the arrow 
shows the direction of E increase. The screening length is A = 30. 

FIG. 7. The radial profiles of local cluster density, p(r), and the energy ratios, E/U{r). The results correspond to the 
different values of attraction energy E = 10,20,30,40,50,60,70,80,90,100 (solid lines), the arrow shows the direction of E 
increase. The screening length is A = 100. 

FIG. 8. The radial profiles of local cluster density, p(r), and the energy ratios, E/U{r). The results correspond to the 
different values of attraction energy E = 10, 20, 30, 40, 50, 60, 70, 80, 90, 100 (solid lines), E = 200(1), 300(2) (dashed lines), the 
arrow shows the direction of E increase. The screening length is A = oo. 

FIG. 9. The scheme of the internal structure of cluster. We display the example of cluster pattern and its local density 
profiles p{r) in radial direction. The intervals of localisation of different zones A,B,C and D are shown, Rn is the compact 
nucleus radius, A is the screening length. 



FIG. 10. The radial profiles of the energy ratios, E/U{r), for the screening length A = 30(a), oo(b). 

FIG. 11. Fractal dimension Df oi cluster branches in the zone D (Fig.^ versus the attraction energy to screening length 
ratio, E/X. 

FIG. 12. The reduced electrostatic interaction energy U* —U/R versus the reduced distance between the centre of compact 
disk- like cluster and unit charge, r"" = r/R ai the screening length A = oo. Here R is the radius of the charged disk. 
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